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ABSTRACT 

Experimental and analytical studies are made of the elastic and 
elastic-plastic response of a three degree of freedom dynamical system 
composed of three rigid and compact masses connected by slender columns, 
the bottom (heaviest) mass being suspended as a ballistic pendulum and 
subjected to an initial velocity step. Such a system is related to the 
behavior of machinery mounted in a ship subjected to an underwater ex- 
plosion. 

The analytical work was performed on a high speed digital computer 
and employed mathematical description of the structural members in which 
force-deflection relations depended upon previous plastic strain history. 

Comparisons are made between experimental and analytical responses 


on the basis of the velocity-time history of the two smaller masses. 


Good agreement is shown. 
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1 Cee rcetuction 


The rresent interest in the resnonse of structures subjected to 
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apicly applied dynamic loads has protuced an increasing amount of 


ryt 


Literature in the area of dynamic loading, The literature covers the 
behavior of structural metals, the respcnse of simple structures, and the 
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theoretical analysis, based esserntialig on energy considerations, was 
able to make good predictions 4° experimental behavior. 

In the case of a similar Multi-story structure, the energy is divided 
between the various roof ani floor masses and it is not possible to use 
energy considerations alone to abtain a theoretical prediction of dynamic 
benavior. Thus, analysis of 2 system above a single degree of freedom 
requires evaluation ot the vaeriois Sguations of notion of the system, Tf 
the structure passes into the plastic region, @ suitable assumption must 
be made to compute the behavior in this region. fhe numerical manipula- 
tions required in solving such a problem become enormous and the use of a 
high speed digital computer becomes necessary. 

This thesis deals with a small two-story structure mounted on a 


large immovable mass, he entire system was reduced to a three degree of 
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freedom system, This particular system can be likened to a piece of 
equipment mounted to a ship's hull. The ship is then subjected to an 
underwater explosion and the system's response is studied, Reference 3 | 
gives examples of equipment and machinery damage and response, wherein 
illustrations of actual shipboard foundation failures are presented along 
with a typical velocity history for a mild shock. The ability to pre- 
dict the response of such a system would aid in the original design of 
such structures. 

A mathematical model was evolved for the above system and a compari- 
son made between the response predicted by the mathematical model and a 
test structure. 

It is felt necessary to mention the fact that the experimental por- 
tion of the project was carried out prior to the evaluation of the mathe- 
matical model and certain complications arose in the mathematical evolu- 
tion that would have been eliminated had the sequence been reversed. The 
values obtained from the mathematical model will be referred to as the 
feeeretical results'' throughout; however, these results cannot be 
truly classified as ''theoretical" since a large amount of experimentation 
with the mathematical model was necessary before a suitable comparison 
could be made. 

The experimental portion cf the thesis was conducted entirely at 
the David Taylor Model Basin, Carderock, Md., during the period 30 June to 
19 July 1961. The theoretical investigation and reduction of data was 
carried out at the U. S. Naval Postgraduate School, Monterey, California, 
during the period 1 January 1962 to 1 April 1962, utilizing the Control 


Data Corporation 1604 computer. 





The experimental investigation will be discussed before the theo- 
retical considerations due to the fact that this was the order in which 
they were actually conducted, A better understanding of the complications 
introduced in the theoretical considerations will be obtained if the factors 


responsible for the complications are presented first. 





2. Experimental Investigation 


Test Structure 





The test structure consisted of two rectangular hot rolled steel bars 
measuring 2'' x 2" x 12", separated and supported by two sets of 1/8" x 2" 
hot rolled steel columns, Figure (la) shows a sketch of th. structure with- 
out the attached measuring equipment. The only variation made to the test 
Structure was the introduction of 3003-H14 aluminum columns for ten (10) 
of the twenty-six (26) test runs. These columns also measured 1/8" x 2", 

The length of the supporting columns, prior to mounting was 12", 
Holes measuring 3/8" in diameter were drilled on the centerline one inch 
from each end of the supporting columns to accommodate a 3/8"-NC bolt for 
clamping the supports to the base and bars mentioned above, 

The rectangular bars and base were drilled and tapped for 3/8''-NC 
threads 3/4" deep at both ends, to receive the securing bolts. Holes 
measuring 3/8" in diameter were also drilled in the sides of the bars for 
attachment of the velocity meters used in the test. Small metal bases, 
fitted to receive the accelerometers used in the test, were welded to the 
upper surfaces of each of the rectangular bars. 

The entire clamping of the supporting columns was accomplished by 
the single 3/8"-NC bolt along with a 1/8" x 2" x 2" clamping plate, see 
figure (lb). The spacing between the base to the center bar and the 
distance between the center bar to the upper bar was approximately 8". 

The single screw attachment of the columns to the base and bars 
proved to be weak as far as giving good fixity at these junctions. This 
poor fixity caused the effective length of the columns to be somewhat in 


excess of the actual free length of 8". Figure (2) is a picture of the 





fully instrumented test structure, the velocity meter coils are on 
the side seen on the picture, The accelerometers are on the upper sur- 
faces facing toward the left. 

The base of the structure contained two 2" x 2" x 2" hot rolled steel 
blocks to be used for attachment of the supporting columns. The base was 
designed to fit a mounting block already attached to the ballistic pendu- 
lum used in the tests, 

Since the width of the supporting columns was large in comparison with 


the thickness, the structure's motion was limited to a single plane. 


Ballistic Pendulum 

The ballistic pendulum was*used as the means for producing the 
impulsive loading on the structure. A picture of the pendulum is shown in 
freure (3). 

The ballistic pendulum consisted of four major components: 1. The 
anvil, upon which the test structure was mounted. 2. The hammer, which 
transferred its energy to the anvil producing the impulsive load, 3. The 
hammer hold-back and release mechanism, 4, The anvil hold-back and 
catching system. 

The hammer consisted of a solid steel cylinder, 8'' in diameter and 48" 
in length. Four smaller cylinders were mounted on the sides to receive 
the supporting wires from the overhead, These cylinders were 4" in da eter 
and 6-1/2" in length. The total weight of the hammer was 785.1 lbs. 

The anvil was likewise 8" in diameter; however, its length was 36". 
The four supporting cylinders measured 3.5'' in diameter and 7" in length. 
The supporting cables were 8 ft in length+1l". The total weight of the 


anvil was 598.6 lbs. 
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ite nawWwer beli-back and vebease mechanism consisted of an electric 
hwist to pull the hammer back co an elevated position, the initial posi- 


scion dependinie unen the impact jesirad, The relessea trrechanism consisted 


i 


> 


of an electro-magnetic letch which could be operated from a remote station 
neay the recording instruments. 

The anvil hold-back and catching system consisted of a ratchet pulley 
and a tie-back line, The line was attached to the envil at one end, passed 


through the pulley and carried » 25 lb. weight et the free end. Upon im- 


Te 


nmact the anvil was forced back and elevate’, Tha 25 1b, weight pulled the 
line through the ratchet pulley and the arvil was prevented from return- 
ing upon reaching its maximum elevetion, 

The ratchet catching system was used to prevent interference between 
the anvil supporting wires and the test structure. The proximity of the 
test structure and the anvil supporting wires, when the anvil was in its 
equilibrium position, was such that the structure came into contact with 
the suprorting wires upon impact. To eliminate this, the anvil was pulled 
back approximately 8-1/2" horizontally prior to impact. 

Since the experimental cesting was conducted prior to a thorough theo- 
retical investigation, the significance of this initial anvil position was 
not realized, and the exact initial horizontal location of the anvil was 
not recorded, Upon examination of the theoretica] aspects it was found that 
this initial position was quite significant and complicated the theoretical 


solution of the problem. More will be said concerning this later, 


Measuring Equipment 
The test structure was instrumented with three bar-magnet velocity 


meters manufactured and calibrated by the David Taylor Model Basin. One 





velocity meter was mounted on each of the three concentrated masses of 

the structure, The meter mounted on the anvil had a sensitivity of 146 
mv/fps, the meter on the center mass (lower’ rectangular bar) had a sensi- 
tivity of 167 mv/fps, and the third meter, mounted on the upper mass (upper 
rectangular bar), had a sensitivity of 151 mv/fps. From tests conducted by 
the David Taylor Model Basin, it was determined that the sensitivities 
varied approximately +3% over a 4 inch stroke. 

Two Statham accelerometers were mounted, one each, on the upper and 
center masses of the test structure. 

The output from the velocity meters went to the Velocity Meter Control 
Unit (Attenuator), manufactured by the David Taylor Model Basin. No ampli- 
fication was possible with this-unit, the maximum output being one times 
the input. This unit also permitted a convenient switching arrangement 
for obtaining a calibration trace on the oscillograph. 

The output from the Control Unit was picked up by a Consolidated 
Electrodynamics Corporation Recording Oscillograph. The oscillograph 
utilized galvanometers of the fluid damped type for recording the velocity 
response on the oscillograph record, The galvanometers had a natural 
frequency of 1000 cps. 

Outputs from’the accelerometers were amplified by Consolidated 
Engineering Corporation 3-KC Carrier Amplifiers. The power supply for 
the above amplifiers consisted of a Consolidated Engineering Corporation 
Oscillator-Power Supply. The amplified accelerometer signals were then 
received by the above mentioned oscillograph. 

A General Radio Company, 1000 cycle, Vacuum-Tube Fork was used to 
produce a time reference trace on the oscillograph record, The oscillo- 


graph also produced a time reference; however, it was felt that the 1000 





eycle trace would provide a geod cheek or tha escillogreph's time scale, 

Since difficulty was encourtered with the interference between the 
test structure and che anvil supporting wires, an extra galwzenometer was 
utilized for recording an indication of contact between the wires and the 
structure. This consisted of a simple battery circuit, using the point of 
probable contact as the switch. When wires and structure made contact, a 
pulse was recorded on the oscillograph record, the duration of this pulse 
gave the period of contact. Figure (4) shows a sketch of the instrumenta- 
tion. 

Appendix VI gives a tabulation of the various instrument serial and 


model numbers, 





Test Procedure 

The following steps were followed for a typical test run 
l. The structure's columns (springs) were placed in position and securely 
fastened, 
2. The anvil was elevated to 2a position where interference between the 
structure and supporting wires would not be encountered. 
3. The velocity meter magnets were adjusted to insure free passage as the 
anvil moved away. 
4, The hammer was elevated to a predetermined position, depending upon 
the step input desired, using the electric hoist. 
5. The recording oscillograph was operated at a paper speed of 4 in./sec. 
and velocity meter and accelerometer calibration traces were recorded, 
6. All calibration settings for the various instruments were recorded. 
7. The instruments were then set for the predetermined settings for the 
actual run, and all values recorded, 
8. The recording oscillograph was operated at a paper speed of 40 in./sec. 
and the hammer released, 
9. The oscillograph was secured when the velocity meter magnets were free 
from the coils or, for the lower input runs, when it was felt sufficient 


data were recorded, 





List otf Symbols 
initial horizontal anvil displacement 
width of rectangular supporting column 
modulus of elasticity 
force produced by flexure of columns (i = 1,2) 
acceleration of gravity (32,2 ft/sec’) 
thickness of rectangular supporting columns 
(subscript) denotes horizontal component 
moment of inertia of the cross section 
spring constant in elastic range 
sprine constant in plastic range 
column length 2 
length of anvil supporting cable 
mass (i = 1,2,3) 
fully plastic bending moment 
axial load on columns produced by structures weight 
tension in hold-back line 
tension in supporting cables 
(subscript) denotes vertical component 
horizontal displacement of M. 


relative spring displacement at which the fully plastic bending 
moment is realized 


vertical displacement of anvil 


dynamic yield strength 
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3. Theoretical Considerations 
Mathematical Model 

The mathematical model of the two-story structure can be visualized 
as a mass - spring system as shown in figure (5a). This idealization 
reduces the structure to a three degree of freedom system. Figure (5b) 
shows the free body diagram for the system. The M, symbols denote the 
masses of the anvil, center mass, and upper mass respectively, the P. 
symbols denote the forces produced by the flexure of the columns, and 


the Xe denote displacements of the respective masses, 


Summing forces for each mass results in the following equations of 


motion, 
x, - — ttn (1) 
= a” (2) 
x, = 4 (3) 


where the force Th is the horizontal component of the force in the 
supporting cables. 

Letting "L" denote the length of the anvil supporting wire, "a" 
the horizontal distance that the anvil is initially displaced, "x'', the 
displacement of the anvil from the initial position "a", and "Y,", the 
vertical displacement of the anvil, we have from figure (6a), 

i = (+x) Sad (2-3) 

where, ae =\Vi2— a@ 


In figure (6b) "W" is the total weight of the anvil plus structure 


(W = 656 lbs.). The symbol "T" denotes the tension in the supporting cables. 





The symbol "'R' denotes the hold-back wire reaction which immediately 


reduces to zero upon impact. From figure (6b) we may write, 


Te = 656 + S88 y, (4) 
Jy —- Sw —- AtX, (5) 
nig a Co6 om ae ne wr 


from figure (6a), 
y= £— ViP-@txy (6) 


( 


O+ x1) x, Bs xe a5 a (8) 
ole a UF; 


L< 
|| 


By substitution of equation (8) into equation (4) and then combining 
with equation (5), we arrive at a value for Th in terms of known or 


computed variables. Substituting this value of Ty, back into equation 


(1) and solving for e produces, 
x +) 
— 1 I u 
X i Sale _ Gt) 656 [1+ ee (9) 


The three equations describing the response of the model are equa- 


tions (2), (3), and (9). 


Basic Considerations 

The test structure followed closely to the portal frame used by 
Kurzenhauser in his tests, the major differences being in the number of 
degrees of freedom and the method of loading. The portal frame had a 
fixed base and the loading tcok place through the upper mass of the struc- 


ture. The two-story structure used in these tests was loaded through its 





base, whith wis attechee to 2 wheihle mass, 

When deforyeticns beyord the elastic lfmit were encountered, it be- 
came necessaty to deviaie fron Wooke's Law and to study the behavior in 
the plastic region. Since the entire system was reduced to a spring - 
mass system it became necessary to assume a force versus deflection relation- 
ship that would satisfy the behavicr of the structure as it passed from the 
elastic region into the plastic region. Figure (7) shows the assumed 
force versus relative deflecticn response. 

The structure's columns behaved elasticaily up to the relative deflec-~ 
tion corresponding to Xu, whereupon the behavior departed from Hooke's Law. 
The simplest way to represent the complicated response in the elastic-plastic 
deformation, is te assum: two linear spring constants, (slope of the lines 
represented by the force deflection curve) one for the elastic range and a 
second for the plastic range. The elastic spring constant is represented 
by k and the plastic spring constant by k*. Upon loading, the columns 
follow the elastic spring constant (bk) until yielding of the columns takes 
place, then they follow the plastic spring constant (k*}. When unloaded, 
the columns follow the fisplaced elastic spring constant (k), forming an 
hysteresis loop as indicated in the figure. 

It is assumed that the transition from the elastic bending moment to 
the fully plastic bending moment in the colimmn cross-sections, occurs in- 
stantaneously at the deflection corresponding to Ku. 

For rectanguler cross-sections, the fully plastic bending moment is 
defined by the relation, 


14 2 
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where G is the yield point or yield strength of the material involved, 


b is the width of the colum, and bh is the thickness of the column, 





A force analysis of one of the columns in the elastic region will 


permit evaluation of the elastic spring constant. Referring to figure 


(8a), and utilizing the cantilever beam deflection equaticns gives the 


following results, 


) (by 
A =- i2/ \2 
= SEL Where: F = total forcing function 
on two supporting 
3 colunns 
ae ame E = modulus of elasticity 
24EI - 


: moment of inertia 


Considering the influence of the compressive force Q, where Q denotes the 
total compressive force on two columns and Ry is the critical Euler column 


loading for two columns, we have, 


_ FL a ca ae ie 
K = jap : . | (10) 


+t? 
or ant 72 rca 


where, 


It has been shown that the approximation used in equation (10, "is 


sufficiently accurate almost up to the critical value of the load, "' [ 6 | 


A force analysis of one of the supporting columns, figure (8b), when 


M, is achieved, gives the following for small deflections 
=Ma = O 
Ze = ae & S 
a Qy + 4+ Mu 





2 a (11) 


Equation (11) is the equation of the straight line representing the 


plastic region shown on the force-deflection diagram of Figure (7). The 
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Spring constant k* for the plastic region is therefore the slope of this 
line, 
= =e (12) 
It will be noted from figure (7. tnat, at the deflectior. Xv the 
spring force represented by equation (10) for the elastic region is equal 
to the force represented by equation (12) for the plastic region. Equat- 


ing these two produces the following, 


Beet wet). Go) aa 
2a ee Me [ B) = - Sx. + aha (13) 


The value for x can be easily obtained from equation (13). 


Complications Introduced by Experiment:1 Set-up 


The two major problems that arose from the experimental model were 
the failure to record the initia! anvil displaeaaiier prior to impact ane 
the lack of fixity of the suppcrting columns. 

Since the initial anvil displacement was not recorded, it became 
necessary to determine this value from the available inforwatior, This was 
accomplished by plotting the sum of the system's total momentum versus time 
for both the exnerimental and theoretical values. The initial displace- , 
ment "a'' used in equation (9) was varied through a reésonable range until 
agreement with the experimental points was obtained. Figure (9) is a 
representative plot. The above mentioned comparison becween theoretical 
and experimental was made for each run where sufficient values of anvil 
velocity were available. Since the length of the velocity trace depended 


upon the time the magnet departed the velocity meter coil, the values for 


15 





the anvil velocity for the plastic runs were not available beyond about 
50 to 100 ™ (millis conds), The value obtained for the initial anvil 
displacement utilizing the above method was consistently 0.7 ft. The 
value arrived at for this displacement should have been fixed for all 
runs, as indicated from the plots, since the anvil was moved to the same 
position for each run during the actual tests. 

The value (a = 0,7') was therefore used for all runs with the ex- 
ception of two of the low input steel runs, which definitely indicated a 
lower value for "a", 

The lack of fixity of the supporting columns introduced a problem as 
to the value for effective length in the formulae for determing k and k*, 
The effective length was varied until phase agreement (peaking at the same 


time) was achieved. The effective lengths found necessary were 9.0 inches 


for steel and 8.75 inches for aluminum, 
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4, Numerical Solution 
Reduction of Equations for Computer 

The first step in reducing the equations to a suitable form for 
introduction into the computer was to reduce the three second order 
equations (2), (3), and (9) to six first order equations. Introducing 
the following dummy variables, 
ae 
Xo 
= OMe 


<x = 4 
t 


then, - 
= X, 


a 


¢ 6@ 


= X, 


CX<e —- C4 - 
iI 


the six first order equations then become, 


Kez (14) 


= 1 (15) 


ee a (16) 


2 


, 25) eis 
a a a ee 656 LI + g ; (17)*- 


M, + (at&)\* , @5e 
2-3.) 3 





s/ = Saar veaian (18) 
2 

-_ a (19) 
3 


* Note that {-JY,, and Y, are given by algebraic relations involving 
X, and f ; see Equations (6) and (7), page 12 
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Letting the symbol D denote the denendent variable, Y denote the 


independent variable, and introducing the following equivalent symbols, 


i = GC = cao 
Y, = X, CK=2£-Y, 
‘ae & YV= Y, 
Y= f 

= % 


Saye AS 
| 


> 


ee 
iH 


Be 


a oS 
i \| 
<0 Ne Cor ox ~S Ch 


Ys 
ex ¥ 


o7 
I 


produces the following equations written in the Fortran language, 


DV) = Y(4+) 
D(2) = Y(5) 
D@) = Y(é) 


(4) = (-F) —(ce/ck) * G56. & (1. + (Y(4) #2 + YV¥"2)/ 
D4) veg EN ine eee / 32.20) *(CC/CK) ¥¥ 2) 


Dts) = (FU) - F(2))/ 0.771 
Die) = Fl2) /0.705 


Appendix I contains the Main Flow Chart, the Force versus Relative De- 


flection Flow Chart, and the actual Fortran program, The program traces 


the time history of the displacements, velocities, and accelerations* 


of the three masses My» Mos and M., starting with the initial conditions 


*Accelerations are given in the actual computer print-out, ho ever, they were 
not used in the analysis and have been trimmed off the sheets shown in Appen- 
dix V. Also, the print statement shown in Appendix I dces not produce dis- 


placements. Displacements were, however, printed in an earlier version and 
it would be very easy to restore them, 
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Numerical Integratior System 
The numerical integration svstem used was the Runge-Rytta method. 
The two major advantages of this system, as stated by Milne [2] » are: 
lL. No special starting procedure is required, 


2. The length of the time interval may be altered during the 
computation, 


The reasens for utilizing the Runge-Kutta method are; fixet of all the 
Fact that the integration relies upon the last computed point, as Listed 
in advantage number one above, :econdly, the routine Was already available 
for immediate use in the U. S. Naval Postgraduate School computer center. 
The system, as with many numerical integration systems, is based on 


® a @ * Ld “ / 
the Taylor Series. For a single first order differential equation y = 


x 
= 


f(t ,y) the following computations are completed to arrive at a value 
y(t + At), 


| = eee” 


K, = Act (4+4t, y+ Kz) 
Y(erat)= 9+ F [K+ 2K, +2Kg + K 
The same basic procedure as above can be applied to a system of first 


order equations such as the six equations shown previously. 
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5. Graphical Cemnerisnn of Tesulte 

This section shews a graphical comparison of the experimental results, 
represented by the circled points, and the theoretical response predicted 
by the mathematical model, represented by the smooth curves, 

The first section displays the results of sixteen runs using steel 
columns, the first seven plots contain a single theoretica} curve since 
these columns experienced no plastic deformation, The next nine plots dis- 
play three theoretical curves corresponding to dynamic yield strengths of 
40,000 psi, 50,000 psi, and 6C,00C psi. The last two plots repeat two of 
the runs extending the time cut-off to 200 MS. 

The second section displays the results of ten runs using aluminum 
columns, The first three runs are in the elastic region and therefore 
display the single theoretical curve. The remaining seven runs display 
two theoretical curves for dynamic yield strengths of 26,000 psi and 30,000 
psi. The theoretical curves shown on pages 44,45 , and 46 actually 
merge, within plotting error, from approximately 90 MS to ICO MS. 

The term "Sten" denotes the initial anvil velocity upon impact (ft/ 


sec). It was Getermined from the velocity meter mounted on the anvil. 
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Section One - Steel 
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6. Discussion of Results 

1. A complete tabulation of the experimental and theoretical points shown 
in the graphical comparison is not included in this thesis. Samples of 

the data are contained in Appendices (II) and (V). The Fortran program used 
is found in Apnvendix (I), and, if desired, may be utilized for obtaining 
additional theoretical results. The experimental results showr on the 


various plots are as accurate as the values used in plotting. 


2, The two maior experimental discrepancies, as noted earlier, were the 
lack of fixity and the failure to record the initial anvil displacement, 
These two items have been suitably discussed in previous sections; however, 
an improvement to the test structure is suggested for future tests. The 
single bolting arrangement used appeared to be inaccquate, as far as pro- 
viding a fixity approaching a cantilever beam, The more rigid connection 
used by Kurzenhauser [1 in his tests consisted of four 5/16" N. C. socket 
head cap screws plus a 1/4" clamping plate. The four screws were located 
near the edges of the clamping plate instead of the centrally located bolt 
used in this thesis, It is therefore felt that the design of the test 
Structure was poor in this respect, and the use of four holding bolts 
would have improved the fixity. 

Nothing further will be said concerning the failure to record the in- 
itial anvil displacement, open that this was vital information that should 
have been recorded. In any future structural tests utilizing the ballistic 
pendulum herein described, the test structure should receive the initial im- 
pact with the anvil hanging vertically (initial displacement zero), thus 


eliminating this initial horizontal force. 


at 











3. The maximum step velocity (initial anvil velocity upon impact) ex- 
perienced by the test structure was 7.83 ft/sec. The ballistic pendulum 
was capable of producing approximately 14.0 ft/sec; however, due to the 
interference between the test structure and the anvil supporting cables, 
the maximum of 7.83 ft/sec. was the highest input before interference was 
experienced, 

4, The limiting value of time used for the graphical comparison of results 
was 100 MS. The reason for cutting off the comparison at this point was 
due to the fact that in reconstructing the initial anvil displacement 

(a = 0.7'), the total momentum vs. time plot was used. This plot was 
limited by the velocity record of the anvil. The velocity meter magnet 

for the anvil departed the coil after approximately 100 MS, up to a step of 
3.86 ft/sec. Above this value for initial velocity, the anvil velocity 
trace was even shorter, averaging about 50 MS, Since the total momentum 
vs. time plot was limited to approximately 100 MS, it was felt that use of 
the predicted anvil displacement beyond 100 MS would not be warrented. 

A second reason for using the 100 MS limit was that the two velocity 
meters mounted on the structural masses produced faulty traces beyond 70 MS 
for the steel columns, at the higher inputs. These faulty traces were lLike- 
wise caused by the velocity meter magnets departing their respective coils. 
The purpose of mounting the accelerometers on the test structure was to 
allow for a continuation of the velocity plots beyond 70 MS by integration of 
the acceleration vs. time traces. It was felt that the increased error due 
to reduction and integration of the acceleration traces would be held to a 


minimum if the number of points added by this method were small. Three 


a2 








points were added, producing a total velocity vs. time trace of 100 MS 
duration. Discrepancies of the order + 5% were noted between velocities 
obtained from the velocity traces directly and by integration of the 
acceleration traces, over intervals where the velocity traces were consider- 


ed valid, 


5. The graphical comparison of results, for steel columns, found on pages 
22 through 37 , consists of apprcximately seven elastic runs and nine 
plastic runs, The experimental points for the elastic runs are compared 
with a single theoretical plot since the theoretical results within this 
range, are not affected by changes in the dynamic yield strength. The plot 
for step = 2,87 ft/sec is a Beidiae ie case. A slight variation between 
theoretical values for = 40,000 psi and CO = 60,000 psi was noted, how- 
ever the variation was quite small (0.3 ft/sec maximum) and therefore only 
the theoretical trace for @ = 40,000 psi is shown. Starting with the value 
step = 3.86 ft/sec, three theoretical curves are shown for the remainder 
of the runs. The curves represent dynamic yield strengths of 40,000 psi, 
50,000 psi, and 60,000 csi. The general trend in the curves indicates the 
increased dynamic yield strength with increased initial velocity. The runs 
for step = 5.02 ft/sec and 7.23 ft/sec do not fit in with the remainder of 
the plots, The dynamic yield strength for these two runs shows a decrease, 
whereas, an increase should be indicated to follow the sequence of values 
previous to and following these respective nee The cause of the dis- 
agreement could possibly be due to an error in the experimental values for 
initial anvil velocity. The comparison indicates a lower velocity input 
than the values used. A change of 0.02" in measurement from the oscillo- 
graph record produces a change of 0.19 ft/sec. Also a variation in support- 


ing columns from one test structure to another could introduce an appreciable 











error, It can be seen in sore of the higher input plots, for steel, that 
the experimental points tend to fall off nore rapidly than would be indicat- 
ed for continuity of plot. As was noved earlier, the last three points 
For the steel plastic runs were computed using the integrated acceleration 
vs, time traces, The added error caused by this procedure may account for 
this discrepancy. 

Pages 38 and 39 show an extension of two of the previous runs out to 
200 MS. The purpose cf these two plots wes to give an indication of the 
agreement beyond the 100 MS limit. As can be seen from the plots, the 


agreement is very good, This further indicates the validity of the computed 


value for initial anvil displacement, 


6. The graphical comparisons using aluminum columns are likewise limited 
to 100 MS. The first three plots show the single theoretical trace whereas 
the remaining seven plets show two theoretical traces, one trace representing 
@G@ = 20,000 psi and a second trace representing @ = 30,000 psi. for 
the three runs represented by step = 2,29, 2.49, and 2.82 ft/sec, the 
theoretical values for @& = 20,000 nsi and & ~ 30,CCO psi merge within 
plotting erret between 90 and 100 MS. 

the experimental values for the velocity-time plot for the aluminum 
columns were taken entirely from the velocity trace since the magnets re- 
mained in the coils somewhat longer than for the steel columns. For the 
two highest input velocities the last three points were checked by integrat- 


ing the acceleration vs. time traces, 


7. Massard and Collins | 2 | , in tests of structural metals under rapid 
loading, demonstrated the time sensitive behavior known as "delayed yield" 


for structural sceel and aluminur, It was found that metals that display 





a definite upper and lower yield point, under slow testing procedures, 
are quite time sensitive and display the delayed yield phenonenon. For 
metals such as aluminum, where straining under slow tests is not dis- 
continuous, the delayed yield effect is not as prominent. From the 
graphical presentation of data, the aluminum plots would tend to substan- 
tiate this conclusion, since the apparent increase in dynamic yield 


strength is somewhat smaller relative to that demonstrated for steel. 


8. An improvement to the testing procedure used would be to replace the 
supporting columns before each run, Obviously, this is a requirement for 
the plastic runs; however, in several instances, the same set of supporting 
columns was used for one very low input elastic run and then a higher 
velocity input plastic run. This procedure only involved about three of 


the runs and was done to save time, 


9. As mentioned earlier, the integration routine utilized in this thesis 
was the RungeKutta method. The time interval used was 1/10 MS. This was 
arrived at by comparison with the results for an exact solution for a 

similar system contained in reference [ 4 | . The only difference between 
the two systems was that the exact solution was based upon a system whose 
base was rigid, The system utilized in this thesis was easily reduced to 


the above system merely by making M, very large. Comparison of the two 


1 
systems in the elastic range (exact solution was valid for elastic response 
only) produced answers that were in agreement to about three significant 
figures. This was considered accurate enough for plotting. 


Since the time interval used in the computer program was 1/10 MS, 


theoretical results could have been obtained for every 1/10 MS; however, 
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this would have produced an excessive amount of output data and the values 
were therefore printed every MS or every tenth time interval. The theo- 
retical curves shown in the graphical presentation of results were plotted 
using points every five milliseconds, 

One inherent difficulty with the Runge-Kutta method, as pointed out 
by Newmark and Chan [5 | , is that “the amplitude will gradually damp it- 
self out even in the case of an undamped system, and the method is there- 
fore not desirable for a long period of time.'' No explanation was given as 
to what constituted "a long period of time''; however, the time limit imposed °* 


in this problem was undoubtedly short enough to prevent this from affecting 


the results. 


10. Perhaps the greatest source of error involved in the experiment was 
the reduction of the experimental data. The data reduction involved first 

fF all a determination of a zero reference level for the various traces. 
This was accomplished by using the edge of the oscillograph record as a 
reference and then comparing this with the 1000 cycle trace for straight- 
ness, In all cases the straight line formed by the 1000 cycle trace agreed 
with the paper edge reference, Amplitudes for the various parameters were 
measured from the respective reference lines using a six inch steel scale. 
A magnifying glass was used for greater accuracy. The scale was read to 
the nearest 0.01". 

As noted earlier, the sensitivity of the velocity meters varied + 3% 
during a 4'' sweep. The response of the accelerometers was good; the only 
drawback in their use was the high frequency pick-up in the first 30 MS. 
This completely obliterated the trace for the first 30 or 40 MS, prevent- 
ing any real comparison of the compatibility between velocity and accelera- 


tion traces, 
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The integratien of the acceleration vs. time plots was accomplished 
by using a K & E 4236 Compensating Polar Planimeter, serial number 6098, 
It is estimated that an error of about 2% was possible using this instru- 
ment since the areas involved were of the order 0.1 sq. in. 

Since the values obtained in the data reduction were accurate to about 
two or three significant figures, the values put into the computer were 
limited ta the same precision, which was also utilized in presenting the 


graphical comparison of results. 


11. {t is to be noted that the parameters chosen for the comparison of the 
results could have been force vs, deflection, acceleration vs. time, or some 
other comparison, The velocity ws. time comparison was chosen because it 
offered the quickest comparison with the least manipulation of numbers, The 
other comparisons mentioned above, would have involved extensive reduction, 


which would have introduced the possibility of increased error. 
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7. Conclusions and Recommendations 

The pratifyingly close corparisons between experimentally observed 
behavior and theoretically predicted behavior of the three mass structure 
considered herein permits the conclusion that the response of similar 
structures may be predicted by mathematical analysis using digital computer 
calculations and involving force-deflection relations for the elastic and 
plastic regimes which invalve prior deformation histery im the ganner des- 
cribed in Appendix I, The analysis in the present case was corplicated by 
a feature of the experimental set-up which would probably have neo counter- 
part in actual practice, There should be no essertial difficulty, in enalys- 
ing cases where the excitaticn cemes from lateral forces having specified 
variations with time rather thea? ee a velocity step of one ot the masses, 
although the latter case is certainly of primary interese tn considering 
response to explosive loadings. ‘Similarly, there should be no o} 
in incorporating further refinevents regarding “large” geometry in the 
analysis should tt be desired to follow the deflection histoiy into scenges 
where the present restricticns to “srall" geometry wight not be sufficiently 
accurate, 

The present experiments Stew that for the higher respsnse velocities 


that accompany the larger input velocity sters, the value of G&G that gives 


Cd 


3 


the best theoretical prediccicr of the response also ircreases. for reasons 
cited by Kurzenhauser [2] »no étteryt was made in this thesic tu make a 
quantitative evaluation of this relation, 

The predicted response Jepends greatly on the details of tha loadc- 
Geflection relations which ceperd, in turn, uron 4a knowledge of the dynamic 
behavior of the material as it is loaded and unleaded in both the elastic 
and plastic regime, [It would be of great advantesge te have bette fauforma- 


tion than seems to be available presently abcut structural materiuls loaded 
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in such a manner; including the variability of @ with strain-rate, 

Future experimental investigations should employ apparatus in which 
a closer control over or knowledge of column fixity is available. Also, 
if a ballistic pendulum type of apparatus is to be used, it is clear that 
the masses should be hanging freely before impact so as to avoid the compli- 
cations introduced in the present experiments by the presence of the hold- 
back cable. Finally, it would clearly be quite advantageous to have a 
single instrumentation capable of reporting the response of the structure 
beyond the 100 MS cut-off without the necessity of augmenting the response 
with a second instrumentation as was necessary in this investigation, where- 
in velocity meter data was extended by use of integrated accelerometer 


data, 2 
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Figure 9 
Plot of total momentum vs. time for an initial velocity input 
(Step) of 1.82 ft/sec. The symbol © denotes experimental 
points and the solid curve represents theoretical values for 
an initial anvil displacement of 0.7 ft. 
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APFENDTX Tf 
FLOW CHARTS AND COMPUTER PROGRAM 
This section contains a very general flow chart for the main 

program, a more detailed flow chart of the force versus deflection 
routine which was the controlling routine of the program, and the actual 
computer program, The force versus relative deflection routine permitted 
the evaluation of the spring forces as the mathematical model experienced 
elastic - plastic response. A list of symbols is included, since several 
symbol changes were necessary in ec: nversion to the Fortran language. The 
figure shown belor and the force-deflection flow chart shown on page 74 


explain the force-deflection subroutine, 


S =e 


—- OCUS OF NEW GENTERS 





The subroutine is er.ered from the main program shown by the box labeled 
MP. Several parameters are picked up such as the relative displacement of 
the columns, denoted by S, the abscissa B of the knee-point, the elastic 
and plastic spring constants k and k*, the value for F where the extension 
of the line represented by k* crosses the vertical or force axis (denoted 


by 4), and the coordinates (P,G) of the present location of the origin of 
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the force-deflection plot. Note that k* is negative as derived on page 15 . 
This is indicated in the figure above. The symbol R is the value of 
displacement used throughout the routine. This value is continually 
corrected for the changing origin. The location of the origin is deter- 
mined from the previous history of the response. The first time that the 
subroutine is entered the values of P and G are zero, that is the origin 
is at (0,0). If the relative displacement S is less than the knee-point 
value B, then the force will be equal to the elastic spring constant k 
times the relative displacement R, where R = S - P or R = S since P = 
zero. When the value of R builds up to exceed the knee-point value B, 

a check is made to see if the column is loading or unloading. The check 
used is a displacement check denoted by the symbol U, where U is equal 

to the present relative displacement S minus the previous relative dis- 
placement SS. As long as the present displacement is larger, the forcing 
function is equa]. to A + k*R. The above routine is followed until the 
columns begin to unload, The present value of relative displacement 
becomes smaller than the previous value and a new origin is computed 

from the relations: P=0O+ H- B and G = k*(P) 

where O is the horizontal value of the previous origin and H is the last 
recorded value of R. The routine loops back to R= S - P where R will 
now be less than the knee-point value B. The value for force will be 
computed from F = G + kR or for values greater than B by F = A + G + k*¥R, 
The above described routine continues until the problem ends and a new 
step input is introduced whereupon all parameters are initialized. The 
same type of routine is also included for the case where the relative dis- 
placement or forcing function or both become negative. This can be seen 


by following the same type of analysis used above. 
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The force-deflection routine is found inside the subroutine labeled 


SUBROU] INE DERIV statement number 77 of the program listing. 
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S518) le 


D 


E = 


CHG - 


S 


S ~ 


TAB - 


H 


List of Symbols 


the initial input velocity to anvil (ft/sec) 
time increment (sec) 
time limit of problem (sec) 
end value of step input (ft/sec) 
initial anvil displacement Taupe (ft) 


correction factor to allow for variable initial anvil 
displacement 


janvil cable length | ee [DD : (et?) 


displacement corresponding to ultimate bending moment 
Eanes 


elastic spring constant , [k| Gian) 
, 
plastic spring constant [k*| (##/in’) 


intersection of plastic force-deflection line with the 
vertical axis (lbs) 


output printing interval (sec) 

step input values for I different values (ft/sec) 
independent variables 

dependent variables 

forcing function in columns (lbs) 

present relative displacement 

previous relative displacement 

time variable 

time expired between print intervals 


last value for relative displacement of origin of force- 
deflection plot 


last value for force origin of force-deflection plot 


*Symbols in | | refer to previously defined terms. 
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SUM 


new force coordinate for force-deflection plot 


new relative displacement coordinate for force-deflection 
plot 


Relative deflection of the columns 
S -—sSS 


total momentum of the structure 
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Main Flow Chart 
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Print - s Time = 


Desired Y Print Interval? 
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Values - Kg 





| Read in new Values 
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#4 Statement numbers are sequential starting with 


PROGRAM KENBOB as statement number one. 
Continued statements count as one statement. 
: | 
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Computer Program for Steel Columns* 
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*The program for aluminum columns is identical except for 


slightly different values for My and Mz. 
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APPENDIX II 
REDUCTION OF DATA 
This section contains a complete set of calculations for one run. 
Also included is a tracing of the oscillograph record for the same run. 
A photograph of the actual record was not possible due to the lack of 
contrast between the trace and the paper. It is to be noted that the 
1000 cycle trace, mentioned earlier, is not shown. This particular 
trace was omitted to prevent cluttering of the other tracings and also 
because of the tedious tracing job required to accurately reproduce this 
trace, 
The data reduction consisted of the following steps: 
1. Establishing reference lines on oscillograph record 
2. Measuring trace amplitudes for velocity parameters 
3. Integrating the acceleration vs, time trace 
4. Computing conversion factors 
5. Converting measured parameters to proper units. 
The run in question was for steel columns and a step input of 5.02 
ft/sec. 
Computation of the necessary spring constants, and other input para- 


meters are also shown, 
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The preceding figure shows traces of the oscillograph record. 
The symbols V denote velocity and A denote acceleration. Direct 


measurements from the figure gives the following values in inches: 


Time v(1) v(2) V(3) 

0 0.88 0 0 

10 * 18 0 

20 * 59 08 
30 * .75 .37 
40 % 75 85 
50 * 86 1.24 
60 * 1.14 1.40 
70 * 1.39 Lae, 


*values not measured 

Calibration pulse measurements: 
V(1l) - 0.720 inches 
V(2) - 0.785 inches 


V(3) 0.820 inches 


A(2) - 1.24 inches 
A(3) - 1.12 inches 
Velocity meter sensitivities and accelerometer calibration constants: 
V(1) - 146 mv/£ps 
V(2) - 167 mv/fps 
V(3) - 151 mv/fps 
**A(2) - 70.7 g 


*kA(3) - 78.0 g 


So 


** values obtained by the David Taylor Model Basin using calibration 
resistance. 
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Instrument settings: 
a. Calibration 


Amplifier gain set at 30 for accelerometers, paper speed 
eet at 4"/sec 


Battery voltage for velocity meters set at 0.6 volts 


b. Run 


Amplifier gain set at 5 for accelerometers, paper speed 
set at 40''/sec 


Combination of the above information gives the following conversion 


factors: 


a. for velocity conversion 


CALIIS. VOLTAGE Cv) ; | 


VrFACTOR SSP REE Eee re — 
. SENSITIVITY CV/fps) PULSE cetT: 


A ce». 
ee ye ee. On 
PFacTOR = 0.146 = CO: 7a in 7 ° in 


° 


ee: = +57 i 


{Pr 
Ve meter = eet i 


b. for acceleration factors 


Ar oo CAL Rs x ONS Vt T RUN GAIN 
PAC hoi = ae : ——————— 
CALT I$ PULSE HIT. CALIB. GAT N 
ile _ 70.7)(22.2) “het 5 gn ty 
Z FACTOR ~~ 1 24- in. rs ZO —_ © SeC=i Nn, 


Az PO me eel 


c. conversion for integration constant 


VER = AN ener OR ss 12. = ractoR ° AREA 
Vac ractoe = 506% y)+ > S25(7¢) . AREA 
Voe- Pe TCk ee = 7.65 © iren Ce) 


a aos _ 9.25 «AREA (Tree ) 
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Areas under acceleration trace: 


2 2 
Time range (MS) Area V(2) (in ) Area V(3) (in ) 
70-80 0.01 -0.04 
70-90 -Q, 14 -0.06 
70-100 Or o2 -0.14 
From the above information: 
Time (MS) V(l) (ft/sec) V(2) (£t/sec) V(3) (ft/sec) 
0 5.02 0 0 
10 * 0.82 0 
20 * Pao AD) 0.39 
30 * 3.43 AS, 
40 % 3.43 4.12 
50 * 3.92 6.00 
60 a walt 6.78 
70 * 6.35 6.63 
80 * 6.43 6.26 
90 * eyes 6.07 
100 * 3690 Doe 


*Not measured 


Evaluation of spring constants - steel 


Column thickness = h = 0.128" 
Column width = b = 1.95" 
Effective length = l ce = 90." 


*modulus of elasticity = E = 29.4 x te psi 
Weight on upper columns = Q, = 23.1 lbs 
Weight on lower columns = Q, = 47.0 lbs 


*Value obtained from actual tension test of sample 
column. See Appendix IIT 
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Evaluation of force-deflection knee points 
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ee Oe eee 
Co ee ae 
Kp = — 5.21 OAK, 
~: an ems) ey, 
KNEE PoinTs 
— Caen OOM ws « 
mee OnGEES a 
Xa. O 6!16 iW 
A = 70O2¢@ Ib, 
(bh) OG = $0,000 pre 
Au, = eres In, 
Au, aes 
A\ = 1906 Ib. 





procedure for 3003-H14 aluminum: 


VW/TE FIX Fe. 


yh. =O 

e) = |.93% 

Neen = TS 
TE 2 to0¥ 10" pse 

Q, - 22 wewslins 

a iy 


*value taken from Alcoa Aluminum Handbook for 3003-H14 Aluminum 


reference | 7 : 
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APPENDIX III 
COLUMN MATERIALS, TENSION TESTS, AND EFFECTIVE MASSES 

Column Materials 

As noted earlier, two different materials were utilized in the test 
colurns; hot rolled steel, with properties similar to those of the steel 
used by Kurzenhauser, and 3003 - H14 aluminum. The aluminum sheet used 
had a yield strength of 21,000 psi and an ultimate strength of 22,000 
psi, as listed by the Alcoa Aluminum Handbook [7] » A 6061 T6 aluminum 
used by Massard and Collins 2 | , would have been more representative 
of the aluminum used in structural work; however, this was not available 
at the time the experimentation was conducted. 
Tension Tests 

Three samples of each of the column materials used were selected and 
shaped into standard tensile test specimens with a 2" gage length to meet 
ASTM Standard E8-54T. A Baldwin Universal Testing machine, along with a 
Riehle extensometer were utilized in making the tension tests. The results 


are tabulated in Table 1 below: 


Table 1 
Steel Aluminum 
Modulus of Elasticity 29.4 x 10° psi 10.27 10° psi 
Yield Point 40,800 psi 20,300 psi 
Ultimate Strength 59,300 psi 22,250 psi 
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Effective Masses 


A very simple approximation of the effective weights of the structure 
was made by including one half of the weight of two supporting columns 
with the upper mass and the weight of two full supporting columns with 
the center mass. The remaining two halves of the lower columns were in- 


cluded with the anvil mass. The figure below shows this procedure. 


eee , 
| MM. PPPECTIVE 
oe. =r oe an x — 
: | 
is . m, LPFECTIVE 
A SE aad i = 
| ay OM, fee re T ne 
a C= ar 
M » 3 | 
ee... — = 


The values of the above effective masses are listed in table 2. 


The units are [pESeen Ee. 


Table 2 
Mass Number Steel effective Al effective 
1 es 18.9 
Zz 0.743 Ov eT Ik 
3 0.717 Os70s 
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APPENDIX IV 


FURTHER PROGRAM MODIFICATIONS CONSIDERED 


Two further modifications to the computer solution were considered. 
One Perl Bieeciion was the determination of the effect of vertical accelera- 
tion upon the structure's columns, The fact that the anvil experiences 
a vertical acceleration upon impact, produces a force in addition to the 
structure's static weight. This vertical acceleration was obtained for 
several runs and found not to exceed about a eylasese This would have 
resulted in a variable correction to the effective value of Q of not 
more than 16%, and the Q/Q, correction was itself rather small, of the 
order of 2%. Thus it was not considered of stfficient importance to 
make this correction which would have greatly complicated ene theoretical 
analysis and would have resulted in considerably more cumbersome dif- 
ferential equations to work with. 

The second modification was an attempt to force the theoretical and 
experimental total momentum vs. time curves to agree exactly, by varying 
the value of initial anvil diepacenent during the run in such a manner 
as to retain the proper noneneen relation. This was tried for several 
runs, resulting in perfect agreement in the plots of total momentum vs, 
time; however, practically no change in the velocity response was noted. 
This modification was likewise not included since the constant value of 


a = 0.7" produced essentially the same results. 
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APPENDIX V 
TYPICAL COMPUTER OUTPUT 
The tabulation of theoretical values for a typical run shows the 
time in MS, velocities V(2) and V(3) in ft/sec, and the total momentum 


in lb-sec designated by the term SUM. 
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APPENDIX VI 


INSTRUMENTATION SERIAL AND MODEL NUMBERS 


Component Mode 1 Serial 
Velocity meter on M, * 159 
Velocity meter on My * Zee 
Velocity meter on M, * 109 
Accelerometer on My A5A- 300-500 4369 
Accelerometer on My A5A- 300-500 4012 
Velocity Meter Control F226 101 
Unit 
Oscillograph 5-119 109025 
Amplifier #1 Wadi lh ays, 198BL10 
Amplifier #2 : 1-113B 641 BD09 
Power Supply 2-105 245R4 
1000 cycle 723-C 2326 


A = = Se 


*None available 
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